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We present a numerical study of persistent currents in quantum rings using current spin density 
functional theory (CSDFT). This formalism allows for a systematic study of the joint effects of both 
spin, interactions and impurities for realistic systems. It is illustrated that CSDFT is suitable for 
describing the physical effects related to Aharonov-Bohm phases by comparing energy spectra of 
impurity-free rings to existing exact diagonalization and experimental results. Further, we examine 
the effects of a symmetry-breaking impurity potential on the density and current characteristics 
of the system and propose that narrowing the confining potential at fixed impurity potential will 
suppress the persistent current in a characteristic way. 
PACS numbers: 73.20.Dx, 73.23.Ra, 71. 15. Mb 



I 

C 

O 

o 



> 

00 

o 
in 
o 
o 
o 

-I— > 



G 

o 
o 



I. INTRODUCTION 

Nanoscopic quantum rings small enough to be in the 
true quantum limit can nowadays be realized experimen- 
tally fj]-§|. Among the quantum effects manifested in 
such systems in the presence of an external magnetic 
field, is the Aharonov-Bohm (AB) effect leading to 
periodic oscillations in the energy spectrum and thus per- 
sistent currents. This phenomenon, first predicted by 
Hund J5|, was discussed in connection with supercon- 
ducting rings [[I)?]] and more recently predicted to occur 
also in one-dimensional metallic rings In the ideal 
case of one electron in a clean, one-dimensional ring, the 
Aharonov-Bohm phase picked up by the electron modi- 
fies the periodic boundary conditions, leading to single 
particle energies given by 



2m \ L ' 



(1) 



where a — 4>/ 4>q is the number of flux quanta penetrat- 
ing the ring and L is the length of the ring. The single- 
particle spectrum is periodic in <f> with periodicity </>o. 
The persistent current associated with state n is 



dE n 



(2) 



However, in realistic systems, interactions, lateral dimen- 
sion, impurities and spin effects complicate the picture. 
In particular, interactions may shift different energy lev- 
els relative to each other, leading to complicated ground 
state patterns with transitions between states with dif- 
ferent spin and/or angular momentum as the Aharonov- 
Bohm flux is increased. In this way, interactions may de- 
crease the period of the oscillations in the ground state 
energy ("fractional Aharonov-Bohm effect"). The first 



systematic study of persistent currents in ideal, one- 
dimensional metallic rings, including temperature- and 
impurity effects but neglecting interactions, was reported 
in [pi. Subsequent approaches include Hubbard model 
calculations [jlO| , use of Hartree- and Hartree-Fock meth- 
ods |llLO|, exact diagonalization studies [OMl5| and 
very recently density-functional calculations jp^,[l7[ . Ex- 
periments in the early nineties reported observations of 
persistent currents in an ensemble of ~ 10 7 Cu rings 
]l8[ , in single gold rin gs |l9| and in a single loop in a 
GaAs heterojunction [BO], all in the mesoscopic range. 
Very recently, Lorke et al. [|1 reported the first spectro- 
scopic data on nanoscopic, self-assembled InGaAs quan- 
tum rings containing only one or two electrons. 

Most of the theoretical approaches mentioned above, 
have the limitation that, e.g. , interactions, spin effects, 
impurities or lateral dimension had to be neglected to 
simplify the calculations. In this paper, we apply the so- 
called current spin density functional theory (CSDFT) 
pl[ including gauge fields in the energy functional. CS- 
DFT was earlier applied to describe the electronic struc- 
ture of quantum dots p2|-p9[. This method, while be- 
ing a mean field approach and thus not exact, has the 
advantage that one can take into account all the above 
effects, being more accurate than Hartree-Fock, and it 
should also be possible to take it to higher particle num- 
bers than the exact diagonalization methods reported 
in the literature. Our aim is, first of all, to exam- 
ine to what extent the CSDFT formalism captures the 
physics due to the Aharonov-Bohm effect, namely the 
periodic variations in the energy spectra as function of 
flux, and the corresponding persistent currents. Thus, af- 
ter introducing the basics of our model in section we 
present, in section [II, the spectra of impurity- free two- 
and four particle rings and discuss how they compare 
qualitatively and quantitatively to corresponding exact 
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diagonalization- and experimental results. In section IV 
we study the effects of a symmetry breaking impurity 
potential on the density profile and persistent current of 
six-electron rings. It is also predicted that, for a fixed 
impurity strength, narrowing the confining potential will 
tend to localize the electrons and suppress the persistent 
current. Finally, section ^ is devoted to discussion and 
conclusions. 



II. MODEL AND NUMERICAL METHOD 

We consider N electrons of effective mass m* , confined 
to a ring with radius Rq by a potential 



V(r) = ~m*u;g (r - R f . 



(3) 



When we examine the system in the presence of an impu- 
rity, we will introduce an additional Gaussian potential 
centered at the bottom of the potential well (x = Rq , y = 
0), 



V/(r) = V exp 



(x - Rq) 2 y 2 
a 2 b 2 



(4) 



The AB flux is provided by a flux tube of constant field 
Bq and with radius r, in the center of the ring; the cor- 
responding vector potential is chosen as 



A r = 



B r/2, r < n 
B r 2 /(2r), r > r t 



(5) 
(6) 



with ri chosen small enough that the electrons themselves 
move essentially in a field-free region. 

Instead of using the quantities N, ui and Ro to de- 
scribe the properties of the system, it is convenient to 
introduce the average II? interparticle spacing r s ,iD (re- 
lated to the average density pip, by r s< iD = V(2p 



ttR /N), and the "degree of one-dimensionality" Cp 
Cp is a dimcnsionless number defined as the ratio be- 
tween the "transverse" (oscillator) gap hut® and the Fermi 
energy; the latter is approximated by the Fermi energy of 
a free, one dimensional Fermi gas with the same density, 
e F = (ir 2 h 2 N 2 )/(8m*L 2 ), where L is the length of the 
ring. Thus, 



Of — 9~z 



(7) 



The higher the value of Cp, the narrower is the ring. We 
will be using values of Cp for which the system is essen- 
tially described by a single channel, but with a smooth 
charge distribution in space. 

For a given set of parameters, we compute the ground 
state charge- and current densities using CSDFT. In this 
formalism, originally introduced by Vignale and Rasolt 



|2l| , one solves the self-consistent Kohn-Sham type equa- 
tions, 



2m* 



2m' 



(p-,4 + A-p)+Vs 



^iS = £i8*i6 (8) 



(We have dropped the arguments r for simplicity). The 
index i labels the eigenstates with spin S = (T,j), and 
A := A + A xc and V s := (e 2 /2m*)A 2 + V s + V H + V xcS 
are the effective vector and scalar potentials. Here, Vh is 
the ordinary Hartree potential and Vs = V+(—)g* iibB /2 
is the external potential, including the Zeeman energy 
(which in our case is set to zero, as the electrons do not 
experience the magnetic field; \xb — eh/ (2m e ) is the Bohr 
magneton). The exchange-correlation vector and scalar 
potentials are 



eA x 



and 



1 f 9 d[pexc(p6,l)) 
p\dy dj 



Vxr.fi = 



d[pe xc {ps,j)} 
dps 



dx 



-i -A 

P 



(9) 



(10) 



where p is the particle density p = p-\ + Pi with ps — 
Si I 2 ■ The paramagnetic current density is given by 
j p - -»V(2m*)Ei 4 [*MV*« - *«V*£], and the real 
current density equals j = j p + (e/m*)Ap. The exchange- 
correlation energy e xc depends on the so-called vorticity 
7 = V x (i p /p)\ z of the wave function. For the de- 
tails of the formalism, we refer to Ref. ]2~]] ], The prac- 
tical computational techniques that we found necessary 
to obtain convergent solutions of the CSDFT mean field 
equations, are given in Ref. |27j| . The parameters and 
results will be given in effective atomic units (a. it.*) with 
energy measured in Ha* — 2Ry* — m* e 4 / (Airheeo) 2 and 
length measured in a* B = h 2 (Anee®) / m* e 2 , where e is the 
dielectric constant and m* the effective mass. The re- 
sults can then be scaled to the actual values for typical 
semiconductor materials. 



III. COMPARISON TO EXACT RESULTS 

In order to check whether CSDFT provides a good de- 
scription of the physics related to Aharonov-Bohm oscil- 
lations, we first apply it to some cases where exact di- 
agonalization results are available, namely impurity-free 
rings containing two and four electrons, respectively. Ex- 
act diagonalization studies of these systems close to the 
ideal, one dimensional limit, were presented by Niemela 
et al. [p^Sjl , who calculated the energy spectra as function 
of the Aharonov-Bohm flux in the presence of interac- 
tions. Starting from the flux-free, non-interacting case, 
one can roughly understand the structure of the full (in- 
teracting) spectra from the following effects: As we have 
seen, the presence of an AB flux induces persistent cur- 
rents and thus favors increasing total angular momentum 
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L (with the corresponding sign) of the electrons. E.g., 
for an even number of particles, as <p is increased from 
zero, L < states come down in energy whereas L > 
are pushed up in energy. The ground state of the non- 
interacting TV-electron ring contains a series of crossings 
between different angular momentum states as the flux 
increases. As interactions are turned on, states with the 
highest possible symmetry in the spin part of the wave 
function are favored due to the gain in exchange energy. 
Thus, for example, triplet states come down in energy 
as compared to singlet states. This may lead to addi- 
tional crossovers, not present in the non-interacting case, 
between states with different spin, in the ground state 
spectrum. Hence, as a first step, we proceed to check to 
what extent CSDFT can produce these features. 
We start by considering a ring with two electrons, choos- 
ing a radius Rq — 1.5 and luq — 1.0. This corresponds 
to the actual values estimated for the experimentally re- 
alized rings recently reported by Lorke et al. [||. Figure 
[l] shows the ground state energy per particle of the two- 
electron system vs. c/>/c/>o as computed from CSDFT. 
(Note that, for this choice of parameters, the electron 
density does not go entirely to zero in the center of the 
ring, so the electrons get partly exposed to the external 
field. This causes the energy spectrum to tilt upwards 
instead of being strictly periodic in 0, in contrast to, e.g. 
, Fig.^, where the ring radius is larger.) 
For the lowest values of tj>, the 5 = 0, L = state is low- 
est in energy, with 5 denoting the spin. Then a crossover 
takes place to the triplet state with L = — 1 and then back 
to 5 = with L — — 2. Exactly the same features were 
found in the exact diagonalization study |l5| , though for 
a more one-dimensional ring. The main difference be- 
tween the two methods is that the cusp corresponding 
to the transition between L = and L = — 2 in the 
5 = state (which is not the ground state here) near 
4>/<po = 1/2 is rounded off in the mean field calculation, 
and the transition between the different L states is grad- 
ual. The reason is the explicit breaking of the rotational 
symmetry in the internal structure of the wave function 
that is mapped out by the self-consistent mean-field so- 
lution. Due to this symmetry breaking, the angular mo- 
mentum is no longer a "good" quantum number, and 
non-integer L-values are allowed. On the other hand, 
cusps at transitions between different spin states are not 
smeared out by the LSDA. 

Furthermore, one can also check by direct comparison 
with the experimental results of Lorke et al. || that even 
quantitatively, our approach provides reasonable predic- 
tions: The energy per particle in the two-electron ring 
at zero external flux can be roughly estimated from the 
experimental data to be around 14meV. Converting our 
calculated result from effective to physical units, using 
the estimated values m* = 0.07m and e = 12.6, the en- 
ergy per particle at zero flux is about 9meV. 
The good agreement between CSDFT and the exact re- 
sults is encouraging and perhaps not too surprising: It 
has previously been shown by Ferconi and Vignale |22] 



that CDFT provides an astonishingly accurate descrip- 
tion of even very small quantum dots, and the same seems 
to be the case here, in the presence of spin. Similar con- 
clusions were reached in a recent, independent, density- 
functional study by Emperador et al. jlffl of the same 
two-electron rings. 
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FIG. 1. Two-particle spectrum: Energy per electron 
vs. flux of the 5 = (solid line) and 5=1 (dashed line) 
states for a two-electron ring with r s> \r> — 2.4, Cf = 18, 
corresponding to the experimental parameters in H . 



Next, we consider a four-electron ring with r Sj iD = 2.5, 
Cp = 10. Fig.|| shows the energy per particle for the low- 
est 5 = and 5=1 states, respectively, as computed 
from CSDFT. We see that the singlet state is always the 
ground state, except for 4'/4 > o near and 1, where the 
two states are nearly degenerate. For both 5 = and 
5 = 1, the total angular momentum changes continu- 
ously from L — at zero flux via L = —2 at 4>/4>o = 1/2 
to L = —4 when one flux quantum penetrates the ring. 
These features agree well with the exact diagonalization 
results in |D| , which again were performed for a more 
one-dimensional ring. 

We note in passing that changing the flux, keeping all 
other parameters fixed, can introduce a weak spin- or 
charge density wave in the ring. 

As before, the main difference between the exact spectra 
and our mean field results is that in our case, the angu- 
lar momentum changes continuously as a function of flux, 
thus smearing out the cusps at transitions between differ- 
ent L-states. This effect, which is due to the approxima- 
tions made in our calculation, is similar to the effect an 
impurity has in the exact spectrum: An impurity poten- 
tial, breaking the rotational symmetry, typically smears 
out the sharp cusps and opens up gaps at level cross- 
ings between different L-states. In the next section, we 
shall study systematically systems with such an explicitly 
symmetry-breaking potential of variable strength. 
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FIG. 2. Four-particle spectrum: Energy per particle 
of the S — (solid line) and S — 1 (dashed line) states 
of a four-electron ring with r s ,iD = 2.5, Cf = 10. The 
numbers along the S — curve indicate the total angular 
momentum L for some values of • 



As an additional test of our method, we have com- 
pared, for given flux, the numerically computed current 
(integrated over a cross-section of the ring) to the deriva- 
tive of the total energy with respect to flux, cfr. Eq.(||). 
We find very good agreement (to three decimal places, 
taking, for example N = 4, r S: io — 2.5, Cf = 10, 
S = — 1 and 4> — 0.2500) with the theoretical prediction 
(g). We thus conclude that CSDFT correctly captures 
the effect of Aharonov- Bohm phases picked up by the 
electrons in the ring. 

For completeness, we also include the spectrum of a 
six-particle ring with r s ^ 1D = 2, Cp = 7 (Fig.||), which 
will be discussed in the next section in the presence of an 
impurity potential. 




IV. EFFECTS OF IMPURITY 



In this section we will study the effects of a single Gaus- 
sian impurity, given by (^), on the persistent current and 
charge density of a six-electron ring. The main feature 
is that the impurity tends to localize the electrons in the 
ring, thus reducing the persistent current. This effect has 
previously been studied by Cheung et al. || in the case 
of an ideal, one dimensional metal ring. We will see that 
our method produces results which are in good qualita- 
tive agreement with those in |9] . 

There is another mechanism which tends to deform the 
electron density: As we shall see, making the ring more 
and more narrow, i.e. increasing the parameter Cf, cre- 
ates a charge-density-wave (CDW) in the ring. In the 
following, we shall study the dependence of the density 
and persistent current on impurity strength and on Cf. 

We start with a direct study of the charge density in a 
ring containing six electrons. Fig Ji] shows the spin down- 
and total densities of a ring with r S) iD — 2, Cf — 7 
and zero flux at various impurity strengths Vo (with 
a = b = 2, sec Eq. (|)). 
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FIG. 4. Spin down- (left column) and total (right 
column) densities of a six electron ring with r S; i£> = 2, 
Cf = 7,4> = 0, S = and impurity strengths (from the 
top) Vo = 0.0, 0.1, 0.5. Increasing Vb gradually localizes 
the electrons and reduces the density at the location of 
the impurity. 



FIG. 3. Six-particle spectrum: Energy per particle of 
the S = (solid line) and S = 1 (dashed line) states of a 
six-electron ring with r s> m — 2, Cp = 7. 



The spin, S, is zero in the ground state, see Fig.||. In 
this case, the density of both spin components is homo- 
geneous and rotationally invariant for Vb = and 0.1, 
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whereas one can see that the electrons start to get lo- 
calized at larger Vq. This localization is accompanied 
by the formation of a spin-density wave (SDW) , with al- 
ternating spin up- and spin down electrons. We also see 
that the density at the impurity center decreases with in- 
creasing Vq. A more systematic analysis of this effect is 
presented in Fig.|5|, where we show the total electron den- 
sity at (x = Rq, y — 0) as function of impurity strength 
for the same ring at three different flux values, and also 
a more narrow ring with Cf = 16, all normalized by the 
density at V = 0. The curves roughly fall on top of 
each other, independent of tfi and Cf and seem to fall off 
exponentially as function of the impurity strength. 




0.2 0.4 0.6 0.8 1 

<£/<J> 

FIG. 5. Density at impurity center vs. impurity 
strength for a set of six-electron rings with r S: io = 2, 
Cf = 7 and <p — 0.05</>o, <p = 0.1</>o and = O.250o, and 
also Cf = 16, <fi = 0.25</>o, all in the 5 = state and 
normalized by their densities at Vq = 0. The shape of 
these curves is roughly independent of flux strength and 
oiC F . 



Fig.|6j shows an example of how the persistent current 
falls off with increasing impurity strength Vq (a = b = 2) 
in a six-electron ring (r s ^\D = 2, Cf = 7) in the S = 1 
state at cf> = O.250o- This choice of parameters, though 
not corresponding to the ground state (see Fig.||), is par- 
ticularly convenient due to the large absolute value of 
the current, making numerical error less significant. Wc 
have checked that other (ground state) sets of parame- 
ters give qualitatively the same behavior. We note that 
the form of this curve, with a plateau at small Vo, fol- 
lowed by steeper fall-off, is very similar to the result by 
Cheung et al. obtained by numerical diagonalization 
of an ideal one-dimensional ring with 20 electrons in the 
tight binding approximation. 
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FIG. 6. Persistent current vs. impurity strength for 
a six-electron ring in the S = 1 state with r s .i£> = 2, 
C F = 7 and <j> = 0.25^0- 



Finally, we examine another effect which tends to local- 
ize the electrons in the ring: It turns out that making the 
ring more one-dimensional, i.e. increasing Cf keeping 
everything else fixed, gradually localizes the electrons, 
creating a strong SDW along with a spatial modulation 
of the total charge density (CDW). This is illustrated in 
Fig.^ where we show density profiles of a six-electron ring 
with r Sl ix> = 2 at zero flux and without an impurity for 
three different values of Cp. 




FIG. 7. Spin down- (left column) and total (right col- 
umn) densities of an impurity-free six electron ring with 
r 8 ,iD = 2, <j> = S = and (from the top) C F = 7, 8, 12. 
Increasing Cf, i.e. narrowing the confining potential, 
gradually localizes the electrons. 
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Such localization and antiferromagnetic ordering of the 
electron spin in a quasi one-dimensional ring confinement 
was recently confirmed by exact diagonalization studies: 
the many-body spectra of quantum rings with up to 6 
electrons could be described by a spin model combined 
with a rigid center-of mass rotation [plf . The stronger 
the ring confinement (i.e. the more narrow the quasi 
one-dimensional ring), or the lower the average particle 
density, the more pronounced are the effects of localiza- 
tion and formation of charge density waves in the internal 
structure of the many-body wave function. 

One might expect that an impurity potential would 
"pin" a charge density wave in the ring, i.e. the persis- 
tent current of a ring with localized electrons should be 
reduced as compared to the non-localized case. We have 
examined the possibility of such a "pinning-depinning 
transition" by computing the persistent current at fixed 
impurity strength as a function of Cp. FigJ| shows the 
result for several different values of the impurity strength 
Vb- We see that the current indeed decreases as the ring 
becomes more narrow; however, there is no "abrupt" 
transition since, as we have seen, localization happens 
gradually with increasing Cf- 
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FIG. 8. Persistent current vs. Cf for a six-electron 
ring with S = 0, r s i£> = 2, <f> = O.250o a t various impu- 
rity strengths, Vb = 0.0, 0.3, 0.6, 0.9 (from the top) and 
a = 2.5, b = 0.5. 



Note the interesting scaling behavior suggested by 
Fig.^: The ratio between any two of the curves is just 
a constant, independent of Cf ; in particular, the ratio 
I(Vq)/I(Vo = 0) is independent of Cf for any Vq. Also 
note that the persistent current decreases with increasing 
Cf even in the zero impurity case, thus making no qual- 
itative distinction between a "clean" and a "dirty" ring. 
This may again be due to the explicit symmetry breaking 
by the LSDA which, as we have discussed previously, in 
a sense mimics disorder even in the case Vb — 0. 



V. CONCLUSION 

The numerical analysis presented here demonstrates 
that the current-spin-density-functional formalism pro- 
vides a suitable tool for describing the effects of 
Aharonov-Bohm phases and impurities in realistic quan- 
tum rings. In particular, we have shown that it repro- 
duces the main qualitative features of the many-body 
spectra and persistent currents, taking into account the 
effects of interactions, spin and deviations from perfect 
one-dimensionality. Furthermore, we have shown that 
the persistent current in the ring may be suppressed 
by narrowing the confining potential at fixed impurity 
strength. 

The main difference between this model and exact 
diagonalization results is that the LSDA introduces a 
breaking of the rotational symmetry in the ring even 
in the impurity-free case. One may hope that this is 
not a problem when describing experimentally realizable 
quantum rings, as a certain amount of disorder and non- 
perfect symmetry is expected to be present in any real- 
istic system. 

With the present experimental progress in fabricating 
and studying few-electron quantum rings, the methods 
described here may turn out useful for suggesting and 
describing future experiments. 
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